arXiv:l 302.02 18vl [astro-ph.CO] 1 Feb 2013 


SUBMITTED TO The As trophy sical Journal 
Preprint typeset using LTpX style emulateapj v. 5/25/10 


DISCOVERY OF A GIANT RADIO HALO IN A NEW PLANCK GALAXY CLUSTER PLCKG1 71 .9-40.7 

Simona Giacintucci 1,2 , Ruta Kale 1 2 3 ’ 4 , Daniel R. Wik 5 ’ 6 , Tiziana Venturi 3 , Maxim Markevitch 5 

submitted to The Astrophysical Journal 

ABSTRACT 

We report the discovery of a giant radio halo in a new, hot, X-ray luminous galaxy cluster recently found by 
Planck, PLCKG1 71 .9-40.7. The radio halo was found using Giant Metrewave Radio Telescope observations at 
235 MHz and 610 MHz, and in the 1 .4 GHz data from a NRAO Very Large Array Sky Survey pointing that we 
have reanalyzed. The diffuse radio emission is coincident with the cluster X-ray emission, has an extent of ~ 1 
Mpc and a radio power of ^ 5 x 10 24 W Hz -1 at 1.4 GHz. Its integrated radio spectrum has a slope of a « 1.8 
between 235 MHz and 1.4 GHz, steeper than that of a typical giant halo. The analysis of the archival XMM- 
Newton X-ray data shows that the cluster is hot (^10 keV) and disturbed, consistent with X-ray selected 
clusters hosting radio halos. This is the first giant radio halo discovered in one of the new clusters found by 
Planck. 

Subject headings : galaxies: clusters: general — galaxies: clusters: individual (PLCKG171.9-40.7) — inter- 
galactic medium — radio continuum: general — X-rays: galaxies: clusters 


1. INTRODUCTION 

Spectacular diffuse radio sources, with very low surface 
brightness and ^Mpc extent, are observed in some of the 
most massive galaxy clusters. Depending on their loca- 
tion, i.e., cluster center versus outskirts, these cluster-wide 
sources are classified as giant radio halos or relics (see 
Ferrari et al. 2008; Venturi 2011; Feretti et al. 2012, for re- 
views; see van Weeren et al. (2012); vanWeeren (2011); 
Kale & Dwarakanath (2010, 2012); Giacintucci et al. (201 1); 
Brown & Rudnick (201 1); Bonafede et al. (2012) for recent 
results). Halos and relics are due to synchrotron emission 
from ultra-relativistic electrons with energy of several GeV, 
spinning in (iG magnetic fields that permeate the hot intra- 
cluster medium (ICM). 

Particle acceleration and magnetic field amplification 
during cluster mergers have been proposed as mecha- 
nisms by which diffuse radio emission is generated in 
clusters. In particular, relics should be produced by 
electron (re)acceleration at merger shocks (Ensslin et al. 
1998; Markevitch et al. 2005; Hoeft & Bruggen 2007; 
Kang & Ryu 2011; Kang et al. 2012), while giant ha- 
los may be caused by reacceleration of lower energy 
relativistic electrons by MHD turbulence during merg- 
ers (Brunetti 2011, and references therein). Although 
alternative possibilities have been proposed for the ori- 
gin of the radio-emitting electrons - the secondary mod- 
els (e.g., Dennison 1980; Blasi & Colafrancesco 1999; 
Pfrommer & EnBlin 2004; Keshet & Loeb 2010; EnBlin et al. 
2011) - turbulent reacceleration models seem to be fa- 
vored by current radio observations (e.g., Venturi et al. 2008; 
Brunetti et al. 2008, 2009; Macario et al. 2010; Donnert et al. 
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2010; Brown & Rudnick 2011; Brunetti et al. 2012). In par- 
ticular, radio halos with very steep spectral index 7 {a > 1.5) 
have been found (the prototype case is A 521 with a 2, 
Brunetti et al. 2008; Dallacasa et al. 2009), consistent with 
turbulent reacceleration in weak mergers, i.e., collisions be- 
tween clusters with relatively low mass (< 10 15 M 0 ) or ac- 
cretion of smaller systems into a massive cluster. Secondary 
models have difficulty explaining such ultra-steep spectrum 
radio halos (USSRHs), due to the large energy required in 
the form of relativistic protons (e.g., Brunetti et al. 2008; 
Macario et al. 2010; van Weeren et al. 2011). However, only 
two giant radio halos are confirmed USSRHs so far - A 521 
and A 697 - based on a well constrained spectrum with five 
data points between 150 MHz and 1.4 GHz (Macario et al. 
2013). Other possible USSRHs still lack confirmation of their 
spectral index (e.g., Giacintucci et al. 2011) and are currently 
under investigation. Therefore, it is necessary to confirm the 
candidates and find more of these USSRHs, with an accu- 
rate determination of their spectral index, before ruling out 
the secondary models. 

A number of dedicated searches for cluster dif- 
fuse radio sources have been carried out in the past 
decade (e.g., Giovannini et al. 1999; Kempner & Sarazin 
2001; Venturi et al. 2007, 2008; Giovannini et al. 2009; 
van Weeren et al. 2009; van Weeren et al. 2011), and the 
number of giant halos and relics has rapidly increased. While 
^50 objects are presently known (including candidates), 
because of their low surface brightness, most of them lack 
detailed spectral information, crucial to test theoretical mod- 
els and understand the origin of the radio-emitting electrons 
(e.g., Venturi et al. 2013). 

The Planck 8 Early Sunyaev-Zel’dovich (ESZ) cluster sam- 
ple (Planck Collaboration et al. 2011b) is an all-sky sam- 
ple of clusters, detected by their multi-frequency signa- 

7 We adopt the convention S u oc i/~ a , where S u is the is the flux density at 
the frequency u. 

8 Planck (http://www.esa.int/Planck), a project of the European Space 
Agency (ESA) with instruments provided by two scientific consortia funded 
by ESA member states (in particular the lead countries France and Italy), 
with contributions from NASA (USA) and telescope reflectors provided by 
a collaboration between ESA and a scientific consortium led and funded by 
Denmark. 


2 


S. Giacintucci et al. 


TABLE 1 

General properties of PLCKESZ G1 7 1.94-40. 65 


RAj 2 ooo (hms) 03 12 57.4 

DECJ 2000 (° ' ") 08 22 10 

Zp e 0.27 

^X[0. 1-2.4] kev (10 45 erg s _1 ) 1.13 

kT (keV) 10.65 

M 50 o (M©) 1.09 x 10 15 

D l (Mpc) 1364.5 

angular scale (kpc/ ") 4.135 


Notes to Table 1 - RAj2ooo and DECj2ooo are the coordinates of the X-ray 
peak, z is from the iron K line. The luminosity Lx [o. 1-2.4] kev> temperature 
kT and mass M 500 are estimated within ^ 500 , where R 500 is the radius cor- 
responding to a density contrast of 500 (Planck Collaboration et al. 2011c). 

ture in the Planck microwave observing bands (30-857 
GHz; Planck Collaboration et al. 2011a). Clusters are de- 
tected through the measurement of the spectral distortion 
of the cosmic microwave background (CMB) due to in- 
verse Compton scattering of CMB photons by the thermal 
electrons in the ICM - the Sunyaev Zel’dovich (SZ) effect 
(Sunyaev & Zeldovich 1972). The Compton y-parameter, a 
measure of the SZ effect, does not suffer dimming due to 
distance, and its integral over the cluster is a measure of 
the total thermal energy in the cluster. Therefore, SZ sur- 
veys are a powerful tool to find massive clusters at high red- 
shifts. A total of 51 new clusters discovered by Planck have 
been confirmed by follow-up X-ray XMM-Newton and op- 
tical observations (Planck Collaboration et al. 2011c, 2012a, 
2012b). Most of these new clusters appear to be hot, mas- 
sive systems with highly irregular and disturbed X-ray mor- 
phologies, and thus good targets to search for new radio halos 
and relics. Recently, double radio relics in one of the new 
Planck clusters, PLCKESZ G287.0 + 32.9, have been discov- 
ered by Bagchi et al. (201 1) using the Giant Metrewave Ra- 
dio Telescope (GMRT). Here, we report a GMRT detection 
of a giant radio halo in the newly discovered Planck cluster 
PLCKESZ G171. 94-40. 65 (hereafter PLCK171). The clus- 
ter, whose general properties are summarized in Table 1, is 
an X-ray luminous (Lx ~ 10 45 erg s -1 ) system at a redshift of 
z Fe = 0.27 =b 0.01, as determined through the XMM Fe K line 
spectroscopy (Planck Collaboration et al. 2011c). An optical 
photometric redshift of z opt = 0.3 1 ± 0.03, based on 29 cluster 
members, is reported by Planck Collaboration et al. (2012a). 

In this paper, we adopt z = 0.27 and the ACDM cosmology 
with H 0 =70 km s -1 Mpc -1 , Q m = 0.3 and Da = 0.7. Under 
these assumptions, 1" corresponds to 4.135 kpc. 

2. RADIO OBSERVATIONS 
2.1. GMR T observations 

PLCK171 was observed with the GMRT at 235 MHz and 
610 MHz in October 2011, as part of a project to search for 
diffuse radio emission in 8 Planck clusters (project 21017). 
The cluster was observed for ^5 hours in dual- frequency 
mode, recording LL polarization at 235 MHz and RR polar- 
ization at 610 MHz. 

The data were collected in spectral-line mode, using the 
GMRT software backend (GSB; Roy et al. 2010) with a to- 
tal observing bandwidth of 32 MHz at 610 MHz divided in 
512 spectral channels. A bandwidth of 6 MHz was used at 


235 MHz. Details on these observations are summarized in 
Table 2, which reports observing date, frequency and total 
bandwidth (columns 1, 2 and 3), total time on source (col- 
umn 4), usable time after data editing (column 5), full-width 
half maximum (FWHM) and position angle (PA) of the full 
array (column 6), rms level (lcr) at full resolution (column 7). 

The data sets were calibrated and reduced using the NRAO 
Astronomical Image Processing System (AIPS) package. The 
data were initially inspected to identify and remove bad chan- 
nels, time intervals and visibilities with radio frequency inter- 
ference (RFI). The data were found to be affected by severe 
phase instabilities caused by ionospheric scintillation during 
most of the observation, leaving only ^ 1 hour of usable time 
at both 610 MHz and 235 MHz. The data were then cali- 
brated. The flux density scale was set using 3C48 and 3C147 
as amplitude calibrators and the Baars et al. (1977) scale. The 
source 0323+055 was used as a phase calibrator. The band- 
pass calibration was carried out using the flux density cali- 
brators. Few of the central channels free of RFI were used 
to normalize the bandpass for each antenna. After the band- 
pass calibration, the central 420 channels were averaged to 35 
channels each 0.78 MHz wide to reduce the size of the data 
set at 6 1 0 MHz, and, at the same time, to minimize the band- 
width smearing effects within the primary beam of the GMRT 
antenna. At 235 MHz, 72 channels were averaged to 18 chan- 
nels of 0.25 MHz width each. 

After further careful editing in the averaged data, a number 
of phase-only self-calibration cycles and imaging were car- 
ried out to reduce residual phase variations and improve the 
quality of the final images. Wide-field imaging was imple- 
mented at each step of the self-calibration process, to account 
for the non-planar nature of the sky. The final images were 
made using the multi-scale CLEAN implemented in IMAGR, 
which results in better imaging of extended sources compared 
to the traditional CLEAN (e.g., Clarke & Ensslin 2006, for a 
detailed discussion see Appendix A in Greisen et al. (2009)). 
We used ^-functions as model components for the unresolved 
features and circular Gaussians for the resolved ones, with in- 
creasing width to progressively highlight the extended emis- 
sion during the clean. The rms sensitivity levels achieved in 
the images at full resolution is ~0. 12 mJy beam -1 at 610 MHz 
and ^0.80 mJy beam -1 at 235 MHz (Table 2). We also pro- 
duced images with lower resolutions, down to 40", by taper- 
ing the uv data using UVTAPER and ROBUST in IMAGR. 
The noise reached in the images at the lowest resolution is 
^0.3 mJy beam -1 at 610 MHz and ~ 2.3 mJy beam -1 at 235 
MHz. Residual amplitude errors should be within ~ 5% at 
610 MHz (e.g., Chandra et al. 2004); a calibration uncertainty 
of ~ 20% is estimated at 235 MHz. 

2.2. VLA observations 

We also reprocessed and analyzed Very Large Array (VLA) 
data at 1.4 GHz from the NVSS (NRAO VLA Sky Survey; 
Condon et al. 1998) pointing containing PLCK171 (project 
AC308). The data were calibrated and reduced in AIPS. We 
used the standard Fourier transform deconvolution method to 
produce the images (CLEAN and RESTORE), and applied 
self-calibration to reduce the effects of residual phase errors 
in the data. The rms sensitivity level achieved in the final im- 
age is ~0.2 mJy beam -1 , with a restoring beam of 58" x 48". 
Our image is ~ 2 times more sensitive than the public NVSS 
image, which has a local rms noise of ^0.4 mJy beam -1 (50" 
restoring beam). Residual amplitude errors are within ~ 5%. 
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TABLE 2 

Details of the GMRT observations. 



Observation 

date 

V 

(MHz) 

Ai/ 

(MHz) 

hot t 

(hours) (hours) 

FWHM, p.a. 
("x",°) 

rms 

(mJy beam -1 ) 

Oct 22, 2011 
Oct 22, 2011 

235 a 
610 a 

6 

32 

5 -1 

5 -1 

15.4x9.0, 64 
5. 5x4. 5, 0 

0.80 

0.12 


Notes to Table 2 - a: observed in simultaneous 235 MHz/610 MHz mode. 


TABLE 3 

Discrete radio sources in PLCK1 7 1 


Source 

RAj2000 
(h, m, s) 

DECj2000 
(0, ") 

5tot, 610MHz 

(mJy) 

Stot , 235 MHz 
(mJy) 

^235 MHz 
^6 10 MHz 

morphology 

SI 

03 12 53.2 

+08 23 12 

11.7 ±0.6 

> 15.5 ±3.3 

>0.3 ±0.2 

NAT 

S2 

03 12 53.9 

+08 23 05 

4.3 ±0.2 

14.6 ±3.1 

1.3 ±0.2 

unresolved 

S3 

03 12 56.9 

+08 22 11 

2.7 ±0.2 

t 

- 

unresolved 

S4 

03 12 57.5 

+08 22 09 

85.8 ±4.3 

t 

- 

NAT 

S5 

03 12 57.3 

+08 21 37 

6.0 ±0.3 

> 8.1 ±1.9 

>0.3 ±0.3 

NAT 

S6 

03 12 58.3 

+08 21 14 

39.2 ±2.0 

66.5 ±13.3 

0.6 ±0.2 

NAT 

S7 

03 12 54.4 

+08 20 35 

1 .9 ± 0. 1 

- 

- 

unresolved 


Notes to Table 3 - Peak coordinates are from the GMRT full-resolution image at 610 MHz (Fig. 1). Flux densities have been measured on the GMRT full- 
resolution images at 610 MHz and 235 MHz (Fig. 2). SI and S5 are partially detected at 235 MHz, therefore lower limits are reported for »Stot, 235 MHz an ^ a - ^ 
Sources S3 and S4 are not resolved in the 235 MHz image. Their combined flux density at 235 MHz is 248.4 mJy. 


3. RADIO GALAXIES IN THE CLUSTER REGION 

Figure 1 presents the GMRT 610 MHz contour image of 
PLCK171 at full resolution (^5"), overlaid on the optical im- 
age from the Digital Sky Survey (DSS). Discrete radio sources 
are labelled from SI to S7. Their flux densities at 610 MHz 
and positions, both measured from the image in Fig. 1, are 
summarized in Table 3. All have optical counterparts on the 
DSS, except for S7. No redshift measurements are currently 
reported for these galaxies in the literature and optical cat- 
alogs. The cluster X-ray centre (Table 1) is approximately 
coincident with the peak of S4. Four sources (SI, S4, S5 and 
S6) show a narrow-angle tail (NAT) morphology. The tails of 
S4 and S6 are both oriented northward, while the tails of SI 
and S5 point toward South-East and North-West, respectively, 
almost opposite to each other. Assuming that the NATs are at 
the redshift of the cluster, the length of their tails ranges, in 
projection, from ^90 kpc (S5) to ~ 1 80 kpc (S6), and their 
radio powers at 610 MHz are of the order of rsj 10 24 -25 W 
Hz -1 . These values are consistent with the range of size and 
radio power typically measured for NAT radio galaxies (e.g., 
Feretti & Venturi 2002). 

The full resolution GMRT image at 235 MHz is presented 
in Fig. 2 (grey scale), with the positions of sources SI to S7 
marked by red crosses and labels. Source S7 is not detected 
at 235 MHz. This is consistent with the fact that the source 
appears very compact at 610 MHz, so it is probably an ac- 
tive AGN with flat or inverted spectrum. The 235 MHz flux 
densities of SI, S2, S5 and S6 are reported in Table 3, along 
with the spectral index a computed between 235 MHz and 
610 MHz. Sources SI and S5 are only partially detected at 
235 MHz, thus the flux density and spectral indeces in Table 
3 should be considered as lower limits. Sources S3 and S4 
appear blended together in Fig. 2, therefore it is not possible 



Fig. 1 . — GMRT full-resolution image at 610 MHz (contours) overlaid on 
the optical DSS image. The restoring beam is 5.5" x 4.5", in p.a. 0°. The 
r.m.s. level in the image plane is 0.12 mJy beam -1 . Contours start at 0.4 mJy 
beam -1 and then scale by a factor of 2. No negative levels corresponding 
to the -3a level are present in this region. Labels indicate the discrete ra- 
dio sources, whose properties are summerized in Table 3. The cluster X-ray 
center is at RA=3hl2m57.4s, DEC=+08d22ml0.3s (Tab. 1). 

to measure their individual flux density at 235 MHz. A total 
of 248.4 mJy is measured for the blend. 
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Right ascension 



0.0010 0.0030 0.0050 0.0070 0.0090 0.0110 


Fig. 2. — GMRT foil-resolution image at 235 MHz (gray scale) with over- 
laid the low-resolution image shown as contours. The restoring beam is 
1 5.4" x 9.0", in p.a. 64° (shown in the bottom-left comer) and 28" x 25", 
in p.a. 73°, respectively. The r.m.s. level is 0.8 mJy beam -1 at full reso- 
lution and 2.3 mJy beam -1 in the low-resolution image. Contour levels are 
6.9, 1 1.5, 18.4, 27.6 mJy beam -1 . No negative contours corresponding to the 
-3 a level are present in the portion of the image shown. The discrete sources 
detected at 610 MHz (Fig. 1; Table 3) are labelled and marked by the crosses. 


TABLE 4 

Properties of the radio halo in PLCK1 7 1 


5235 MHz (mJy) 

483 ±110 

5*61 o MHz (mJy) 

> 50 

5*1400 MHz (mJy) 

18 ±2 

a (235MHz -A 1400MHz) 

1.84±0.14 

A 400 MHz (W Hz -1 ) 

(4.9 ±0.1) x 10 24 

Linear size (Mpc) 

- 1 


A low-resolution (28" x 25") image at 235 MHz of 
PLCK171 is presented as contours in Fig. 2, overlaid on the 
full-resolution image in gray scale. The image shows the 
residual emission after subtraction of the CLEAN compo- 
nents of sources SI to S6 from the uv data. The CLEAN com- 
ponents were obtained from an image produced using only the 
baselines longer than 3 kA, containing information on struc- 
tures on angular scales <1.4'. The presence of a large diffuse 
source filling the central region of the cluster is evident. Its 
large linear size 800 kpc in Fig. 2), roundish morphology 
and the central cluster location are consistent with the proper- 
ties of a giant radio halo. 

We note that the peaks of the radio halo image in Fig. 2 
are approximately coincident with the position of some indi- 
vidual sources. To estimate the possible residual contribution 
of these sources, we made an image of the halo region from 


the source-subtracted data set, using only the long baselines 
(> 3 kA). We indeed found some residual emission possibly 
associated with the blend of S3 and S4 (~ 7 mJy) and with S6 
6 mJy). This will be taken into account when measuring 
the flux density and spectral index of the halo in Sect. 5. 

Hints of the halo are already visible on the public NVSS im- 
age at 1 .4 GHz. To obtain a better quality image of the diffuse 
source, we re-analyzed the NVSS pointing (see Sect. 2.2). 
Our 1.4 GHz images are presented in Fig. 3. Panel a shows 
the total radio emission, i.e. the halo and the radio galaxies. 
Following the procedure adopted at 235 MHz, we re-derived 
the image using only the baselines longer than 0.8 kA (angular 
scales < 4') to identify the contribution of the discrete radio 
sources in the cluster area. A central region of the resulting 
image is shown in panel b, with the position of the sources 
detected at 610 MHz (Fig. 1) highlighted by crosses and la- 
bels. Even though the resolution of the two images is very 
different (~ 5" at 610 MHz and ~ 50" at 1.4 GHz), there is a 
good match between the location of the radio sources at 610 
MHz and the structure imaged at 1 .4 GHz. The CLEAN com- 
ponents associated with this image were then subtracted from 
the uv- data, and the resulting data set was used to obtain an 
image of just the radio halo, shown in Figure 3c (colors and 
contours). An overlay of this 1.4 GHz diffuse emission on 
the GMRT 235 MHz image of the radio halo, convolved with 
a circular beam of 40", is shown in Figure 3d. The overall 
shape of the halo is similar at both frequencies, while its ex- 
tent is slightly larger at 1.4 GHz (~ 1 Mpc). The difference is 
likely due to the low sensitivity and quality of the 235 MHz 
observation, which was strongly affected by ionospheric scin- 
tillation during most of the observing time (Sect. 2.1). 

5. INTEGRATED SPECTRUM OF THE RADIO HALO 

Using the images of the diffuse component alone (Fig. 2 and 
Fig. 3c), we measured a total flux density of the halo of 483 d= 
110 mJy at 235 MHz 9 and 1 8 ± 2 mJy at 1 .4 GHz in the region 
defined by the first contour in Fig. 3d (3cr). These are also 
reported in Table 4, along with the other main observational 
properties of the halo, and in Fig. 4. The derived total radio 
power at 1.4 GHz is P\agrz = (4.9 ± 0.1) x 10 24 W Hz -1 . 

Our flux density values at 235 MHz and 1.4 GHz corre- 
spond to a very steep spectral index, a= 1.84 ± 0.14, which 
makes PLCK171 one of the steepest-spectrum giant halos 
known so far. However, large uncertainties affect our mea- 
surement. In particular, it is possible that some of the halo 
flux density is missed in both the reprocessed NVSS image 
and the GMRT 235 MHz image due to the short duration of 
the observations, which results in sparse coverages of the uv 
plane and loss of sensitivity to structure. Furthermore, the low 
resolution of the 1.4 GHz image does not allow an accurate 
removal of the discrete radio galaxies in the halo region, and, 
consequently, the halo flux density could be overestimated. 

A flux density of ^45-50 mJy is measured from low- 
resolution (~ 30-40") images at 610 MHz (not shown here), 
obtained after the subtraction of the discrete sources. How- 
ever, these images detect only the brightest peaks of the halo 
and, therefore, its flux density at this frequency is only a lower 
limit. The 610 MHz measurement is ~ 40% lower than the 
value expected based on the spectral index between 235 MHz 
and 1.4 GHz (Fig. 4). Our experience with the GMRT at 610 
MHz suggests that the flux density losses at this frequency 

9 The residual 13 mJy, possibly associated with the individual sources (see 
Sect. 4), have been added into the error. 
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0.00018 0.00054 0.00090 0.00126 0.00162 0.00198 0.0010 0.0046 0.0082 0.0118 0.0155 0.0191 

Fig. 3. — VLA-D array images at 1.4 GHz of PLCK171 obtained from the re-processed NVSS data. Contour levels are -1 (dashed), 1,2,4, 8, 16, ... x 3cr. a: 
Contours and color scale image of PLCK171. The restoring beam is 58" x 48", in p.a. -5° and 1 <7=0. 17 mJy beam -1 , b: Image obtained using only baselines 
longer than 0.8 kA. The restoring beam is 50" x 42", in p.a. -7° and 1 <j= 0.4 mJy beam -1 . Crosses and blue labels show the position of the discrete radio sources 
detected at 610 MHz (Fig. 1, Table 3). c: Image of the diffuse radio halo after the subtraction of the discrete radio sources derived from the image in panel b. 
The image has been restored with a circular beam of 50". The r.m.s. noise is 1 <t = 0.12 mJy beam -1 . d: 1.4 GHz contours (same as panel c) overlaid on the 235 
MHz image of the radio halo convolved with a circular beam of 40". 


can indeed be conspicuous, up to ~40% (Venturi et al. 2008; 
Dallacasa et al. 2009; Macario et al. 2010). 

6. X-RAY ANALYSIS 

PLCK171 was observed with XMM-Newton for 14 ks in 
2010 August under ObsID 0656201101. The EPIC data 
were reduced using the Extended Source Analysis Software 
(XMM-ESAS) 10 package part of SAS version 1 1 .0.0. For the 
EPIC MOS detectors, the analysis of galaxy clusters with this 
software was introduced by Snowden et al. (2008), but the 
methodology has since been expanded to include the EPIC 
pn data (e.g., Bulbul et al. 2012, for a more detailed descrip- 
tion). In summary, a filtered event list is created by removing 
periods of excessive count rate in the light curve caused by 
proton flaring events (retaining 13.7 ks, 13.4 ks, and 8.2 ks 
of the exposures for the MOS1, MOS2, and pn data, respec- 

1 0 http ://heasarc . gsfc . nasa . gov/ docs/ xmm/ xmmhp xmmesas . html 


tively), excluding anomalous MOS CCDs (in this case CCDs 
4, 5, and 6 for MOS1), and identifying and excluding point 
sources. Spectra and images are then extracted, including for 
the quiescent particle background, which is derived from a 
database of filter wheel closed observations that are matched 
to data from the unexposed corners of the chips. The remain- 
ing background components (residual soft proton contamina- 
tion, instrumental lines, solar wind charge exchange (SWCX) 
emission, and the cosmic fore- and backgrounds from the lo- 
cal hot bubble, the Galaxy, and extragalactic unresolved point 
sources) are all explicitly modeled and determined empiri- 
cally when fitting spectra in XSPEC. While robust best fits 
of these components can be trickier to obtain in shorter expo- 
sures such as this one, the smaller angular extent of this cluster 
provides for a large local background region from which their 
parameters can be sufficiently constrained. 

In this case, nearly all of the cluster emission is taken to 
fall within a circle of 320" radius; EPIC MOS and pn spectra 
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Fig. 4. — Spectrum of the radio halo between 235 MHz and 1.4 GHz. 


from this region and the remaining field of view out to 800" 
(containing only background) are simultaneously fit. We find 
a global cluster temperature of 10. 1 db 0.7 keV, in good agree- 
ment with that reported in Planck Collaboration et al. (201 lc) 
of 10.65 ±0.42 keV. Our slightly lower temperature likely 
reflects a different source region and background treatment, 
and our greater uncertainty is a consequence of treating back- 
ground components as free parameters. For the image analy- 
sis, the best-fit values for the non-cosmic backgrounds (qui- 
escent particle, soft proton, and SWCX) are translated into 
images and subtracted from the data. 

The resulting adaptively smoothed, full band (0.4-7. 2 keV) 
image is presented in Figure 5a with the 1.4 GHz radio halo 
contours overlaid (same as Fig. 3c). Structurally, the central 
region is elongated in the North-West/S outh-East direction 
with the X-ray peak slightly towards the SE side from the 
large-scale emission centroid. The radio halo and X-ray sur- 
face brightness generally follow each other, as often observed 
in clusters hosting radio halos. The asymmetric ICM mor- 
phology is consistent with a recent merger event. 

To investigate the merger state in more detail, we con- 
structed a temperature map using smoothed narrow-band im- 
ages fitted with the MeKaL plasma model, using a method de- 
scribed in Markevitch et al. (2000). Images are extracted in 
seven bands (0.4-0.8 keV, 0.8- 1.1 keV, 1.1-1.35 keV, 1. 9-2.1 
keV, 2.35-3.2 keV, 3 .2-4.8 keV, and 4.8-7.2 keV) that avoid in- 
strumental lines. A local background, accounting for Galac- 
tic and extragalactic non-cluster sources, is estimated from 
the region outside the cluster in each image and subtracted, 
after which the image is corrected by vignetting and expo- 
sure maps. Using the ESAS task comb, which adjusts each 
exposure map to the response of the MOS2, medium filter re- 
sponse, the three EPIC instruments are combined to maximize 
the signal-to-noise before smoothing. The images, with 2.5" 
pixels, are variably smoothed by a 1 0 pixel Gaussian kernel 
at the image edges, which transitions to a ~ 5 pixel kernel in- 
versely following the wide band surface brightness to the 0.2 
power. Chip and point source gaps are interpolated over (for 


purely cosmetic purposes). Then, an absorbed MeKaL model 
with fixed absorption (2.8 x 10 21 cm -2 ), abundance (0.3 so- 
lar), and redshift (0.27) is fit to the coarse spectrum in the 
seven energy bins corresponding to each pixel. The resulting 
best-fit temperatures are shown in Figure 5b, which have 1- 
sigma errors on the order of 1 keV in the center and ~ 2 keV 
near the edge. 

The general picture that emerges is a roughly isothermal 
central region, without a central cool core and with a slight, 
but only marginally significant, temperature peak (^11 keV) 
coincident with the X-ray peak. The temperature peak extends 
in the same direction, but more narrowly, as does the broad 
band surface brightness. This elongation follows that evident 
in the core of the radio halo, suggesting a merger along the 
NW-SE axis. The temperature of the cooler region to the SW 
of the X-ray peak (~ 7 - 8 keV) deviates from the tempera- 
tures along the NW-SE axis at the 2.3 a level and is confirmed 
in direct spectral fits. Deeper X-ray observations are neces- 
sary to discern exactly what this feature is and its relationship 
to the merger. 

7. DISCUSSION 

The Planck mission has discovered several new clusters us- 
ing the SZ signal. These systems were not detected from the 
all sky X-ray surveys, such as the ROSAT All Sky Survey, be- 
cause of their low surface brightness and being near the flux 
limit of the surveys. The newly discovered Planck clusters are 
mostly massive systems with highly disturbed ICM. With the 
current knowledge of the connection between cluster merg- 
ers and the presence of giant radio halos (Cassano et al. 2010, 
and references therein; see also Rossetti et al. (2011)), these 
new clusters are excellent targets for searching for new halos. 
In this paper, we reported the discovery of the first giant ra- 
dio halo in these new Planck clusters - in PLCK171 - based 
on GMRT observations at 235 MHz and 610 MHz and re- 
analysis of the VLA 1.4 GHz NVSS data. 

The diffuse halo is ~1 Mpc in extent and is cospatial with 
the brightest X-ray emission, as typically observed in radio- 
halo clusters (e.g., Govoni et al. 2004; Venturi et al. 2013). 
The spectral index of the halo between 235 MHz and 1400 
MHz is quite steep, a=1.84±0.14. This is one of the steep- 
est slopes measured for a giant radio halo so far - in fact, 
it is similar to a ~ 1.9 of the prototype USSRH in A 521 
(Dallacasa et al. 2009). However, the spectral index of the 
PLCK171 halo is affected by large uncertainties. In particu- 
lar, due to the short usable time of our 235 MHz observation 
and short integration time of the NVSS pointing, it is possible 
that the halo extends more than imaged here, and a fraction 
of its flux density is missed at both frequencies. Moreover, 
the low resolution of the 1.4 GHz image does not allow an 
accurate subtraction of the discrete radio galaxies enclosed in 
the halo emission. Therefore, deeper observations at 1 .4 GHz 
with higher angular resolution, and deeper exposures at lower 
frequencies are essential before any conclusive interpretation 
of the spectral index of this halo can be made. If future ob- 
servations confirm its steep radio spectrum, PLCK171 would 
be another important case of USSRH (of which only a few 
other examples are known so far; e.g., Brunetti et al. 2008, 
Macario et al. 2010, van Weeren et al. 2011, and a candidate 
in Giacintucci et al. 2011), whose study is crucial to under- 
stand the physics behind the origin of radio halos in general. 
While such objects are expected in turbulent reacceleration 
models, the existence of halos with a > 1.5 poses serious 
problems for secondary models due to the required large en- 
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Fig. 5. — a: Adaptively smoothed XMM-Newton image in the 0.4-7.2 keV band, b: Temperature map. The colorbar displays the temperature in keV. In both 
panels the 1.4 GHz image of the halo is overlaid as contours (same as Fig. 3c). 
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Fig. 7. — Distribution of the radio-halo clusters in the A .4 GHz - Az(< 
57?5 oo) sr Mpc 2 plane (Basu 2012). The green symbol marks the position 
of PLCK171. The red symbols are the known radio-halo clusters, the yellow 
symbol is the mini-halo in A 2390 and blue symbols are the upper limits from 
the GMRT radio-halo survey (Brunetti et al. 2007). 


Fig. 6. — Distribution of the radio-halo clusters in the P\aghz~Lx plane 
from Brunetti et al. (2009). Filled black symbols are the GMRT radio-halo 
clusters and open black symbols are other radio-halo clusters from the liter- 
ature. The magenta point marks the position of PLCK171, using the cluster 
X-ray luminosity estimated wit hin R 5 qo (Planck Collaboration et al. 2011c) 
and the 1.4 GHz radio power in Table 4. The solid line is the best fit to the 
distribution of giant radio halos from Brunetti et al. (2009). 

ergy in cosmic ray protons (Brunetti et al. 2008, Macario et 
al. 2010, van Weeren et al. 2011). Such energies are above 
the current upper limits from the 7 -ray observations of nearby 
clusters (Aharonian et al. 2009, Ackermann et al. 2010, Jel- 
tema & Profumo 2011), assuming the fiG cluster magnetic 
field values indicated by Farady rotation measure studies (e.g., 
Bonafede et al. 2010 and references therein). 

The analysis of the X-ray XMM-Newton observation of 
PLCK171 reveals high temperature (^10 keV) and signa- 


tures of a recent merger in both the surface brightness and 
gas temperature distributions. A NW-SE elongation of the 
central X-ray emission, coupled with a similar structure in 
the temperature map, suggests that the merger occurred along 
this axis. The radio halo displays a similar asymmetry in its 
central region, suggesting a correlation between the thermal 
and non-thermal emissions, as typically observed in radio- 
halo clusters. 

Figure 6 shows the distribution of the clusters of the GMRT 
radio-halo survey (Venturi et al. 2007, 2008), along with other 
radio-halo clusters from the literature, in the P\agrz-Lx plane 
(from Brunetti et al. 2009). The bimodality evident in this 
diagram (with some clusters lacking diffuse radio emission, 
while those that exhibit radio halos follow a correlation) is 
one of the main outcomes of the GMRT radio-halo survey and 
provides quantitative support to reacceleration models for the 
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Fig. 8. — Spectral index of radio halos in the 235 MHz-1.4 GHz interval as 
function of the cluster X-ray temperature (adapted from Venturi et al. 2013). 
The magenta points marks the position of PLCK171. 

origin of radio halos (Brunetti et al. 2007, 2009). The separa- 
tion between clusters hosting a radio halo and those without 
(upper limits) is most likely caused by the different dynamical 
properties of the clusters: clusters without a radio halo tend to 
be more relaxed systems, while radio halos are exclusively 
found in disturbed clusters (Cassano et al. 2010). 

PLCK171 (magenta symbol) falls on the radio halo corre- 
lation. Furthermore, our X-ray analysis shows that the cluster 
is unrelaxed, as expected for a system hosting a radio halo. 
Ongoing dynamical activity in PLCK171 is also suggested 
by the presence of four NATs in its central Mpc region, with 
tails pointing away from the cluster center, suggesting infall 
of multiple subctructures onto the main cluster and bulk mo- 
tions of the ICM driven by the ongoing merger (e.g., Bliton et 
al. 1998). Thus, the radio halo in PLCK171 provides further 
support to the existence of a direct link between the radio halo 
phenomenon and cluster mergers. 

A scaling relation between the radio power at 1 .4 GHz and 
the integrated SZ effect measurement (Y) has been recently 
presented by Basu (2012). Such a relation is a direct probe of 
the connection between the mass of the cluster and presence 
of a radio halo. In Fig. 7 we report the position PLCK171 
in the P\aghz-Y plane using the integrated Y published in 
Planck Collaboration et al. (2011c). PLCK171 (green square) 
follows the scaling relation obtained for the other radio-halo 
clusters. 

What is unusual about PLCK171 is its combination of high 
gas temperature (~ 10 keV) and ultra-steep spectrum of its 


radio halo (a ~ 1.8, with the caveats given above). If the steep 
spectrum is confirmed, it would be the hottest cluster with a 
USSRH - all other USSRHs have kT ~ 5-9 keV, as visible in 
Fig. 8, that reports the distribution of radio-halo clusters in the 
a 1 400 MHz ~~ kT plane adapted from Venturi et al. (2013). The 
magenta point marks the position of PLCK171 using kT = 
10.1 ± 0.7 keV, as measured in Sect. 6. Thus, this object may 
provide an important data point for our understanding of giant 
radio halos in clusters. 

8. CONCLUSIONS 

We have discovered a ~ 1 Mpc-size radio halo in 
PLCK171, using GMRT observations at 235 MHz and 610 
MHz and re-analysis of NVSS 1.4 GHz data. This is the 
first giant radio halo found in the new clusters discovered by 
Planck. 

With a spectral index a ~ 1 .8, this source might be another 
system with a very steep radio spectrum, so far found in only 
a few other clusters. However, deeper multi-frequency obser- 
vations are required before any definitive interpretation of the 
spectral index of this radio halo. 

The X-ray luminosity of the cluster and radio power of the 
halo at 1 .4 GHz are consistent with the correlation known to 
be followed by radio-halo clusters. Its location in the P\agsz- 
Y plane also follows the scaling relation found for clusters 
with radio halos. The asymmetries observed in the X-ray sur- 
face brightness and temperaure distributions, derived from the 
XMM data, indicate a disturbed ICM resulting from a recent 
merger. This is in line with the expectation that giant radio 
halos are exclusively associated with cluster mergers. 

With an X-ray luminosity of 1.13 x 10 45 erg s -1 and a red- 
shift of z = 0.27, PLCK171 represents the kind of clusters that 
have been missed in X-ray flux limited surveys, such as the 
Rosat All Sky Survey, and, consequently, in radio surveys 
based on X-ray selection criteria, such as the GMRT radio- 
halo survey (Venturi et al. 2007, 2008). It is therefore im- 
portant to consider biases due to clusters such as PLCK171 in 
X-ray flux limited samples selected for radio surveys. 
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